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Abstract 
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1 Introduction 


It is becoming increasingly clear that the study of strongly coupled gauge theories in unphysical 
settings, such as finite volume or fixed gauge field topology, can be extremely fruitful. Prom the point 
of view of lattice simulations, the restriction to finite volume is of course ideal. It was realized long 
ago that it can be advantageous to consider the propagation of the pseudo-Goldstone bosons inside 
volumes too small to even contain one Compton wavelength 1 /m, of these light degrees of freedom 
Si- This idea was explored in depth in the chiral Lagrangian framework of full QCD ||, Q, as well 
as in a class of sigma models j|] . Already at that time numerical simulations successfully reproduced 
these predictions for sigma models [j6l, j?|. 

In the meantime lattice simulations of QCD with very light quarks are beginning to become feasible. 
If one restricts oneself to the same regime of finite size L = H 1 / 4 , where m n <C 1 /L, it was noted by 
Leutwyler and Smilga that gauge field topology now plays a much more important role than in the 
infinite-volume theory §]. This fact can be an advantage also in numerical simulations. If one has 
very definite and different predictions in the various sectors of fixed topological charge v, it means 
that one can perform a whole series of independent fits to lattice data, and not just one. 

It is thus of interest to combine finite volume with fixed gauge field topology, and to consider QCD 
correlation functions in such a situation. A systematic framework for analysing QCD in the regime 
m n « 1/L < 47tF, the so-called e-expansion 0, turns out to generalize easily to sectors of fixed 
topology. In this paper we shall describe our computation of finite-volume meson correlation functions 
in the quenched analogue of this e-expansion. We shall also provide new results for the unquenched 
theory in the same finite-volume regime, but restricted to fixed topology. The completely quenched 
theory is well-known to be special. Eliminating the fermion determinant renders the theory non¬ 
unitary. Still, it is obviously of great interest to understand this truncation of the real theory in 
view of lattice simulations. Here we shall adopt the point of view that the quenched theory can be 
given some meaning close to the full theory, over a certain range of scales. The difficulties of the 
finite-volume quenched theory in the chiral limit at the quark level show up again in the effective field 
theory description, and there is no way around this fact. 

Our paper is organized as follows. In section ^ we discuss the two different ways in which to do 
quenched chiral perturbation theory, the supersymmetric method, and the replica method. Section 
discusses two different regimes of chiral perturbation theory at finite volume, yielding two different 
perturbation series known as the p- and e-expansions. We point out the peculiarities of the counting 
of orders in these two expansions when going to the quenched theory, and explain why sectors of fixed 
gauge held topology can be particularly useful in the quenched case. In section f| we recalculate the 
leading order correction to the quenched chiral condensate in the e-expansion, and show that it gives 
the same result as computed earlier in the replica formalism in ||. Section |E| contains the results 
of our quenched calculation of mesonic correlation functions at fixed topology, to leading order in 
the e-expansion. We illustrate our results for typical values of lattice parameters, and point out how 
these results can be used to determine the low-energy parameters of the chiral Lagrangian such as the 
pion mass m n and the pion decay constant F. We also present various checks on our calculations, 
such as consistency with exact Ward identities. In section ^ we present the results of the analogous 
calculations for full QCD with Nf light quark flavours at fixed gauge field topology. Section ^ contains 
our conclusions, and in two appendices we collect various technical details relevant to our calculations. 
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2 Supersymmetry and Replicas 


There are two different, yet equivalent, methods for performing quenched chiral perturbation theory. 
The standard one is based on a chiral Lagrangian extended to a supergroup [jil]], and the other on 
the replica method [JlTJ . In the former, internal quark loops are cancelled by a mechanism similar 
to the cancellation of vacuum diagrams in space-time supersymmetric field theories. In the latter, 
quark loops are removed by setting the number of quarks equal to zero. The equivalence of the 
two methods in chiral perturbation theory follows quite easily once one sees the two different sets of 
Feynman rules. In both cases quenching simply corresponds to removing the fermion determinant in 
the original theory. 


Let us begin at the quark level. Because we shall be interested in computing quantities that include 
quarks on external lines, we denote by N v the number of such “valence quarks”. In the supersymmetric 
formulation we include appropriate sources J for these valence quarks (for simplicity restricted to quark 
bilinears here), but cancel their contributions, at vanishing sources, by quarks of wrong statistics. The 
generating functional is thus 


Zsusy[J] — / [^4/J 


det (iL/> - m v + J) _ S [ A] 
det (0 — m v ) Nv 


( 2 . 1 ) 


where the determinant in the numerator is over N v fields. The resulting theory is not space-time su¬ 
persymmetric, because there are no superpartners of the gluons, but internal quark loops are cancelled 
by their ghost partners. 


In the replica method we embed the N v valence quarks in a theory with N quarks in total (all of 
which are of ordinary fermionic statistics). At that stage the generating functional reads 


ZRepiica[J] = JidAfj] det(0-m v + J) det(i$ - m v ) N Nv e 5[A] 


( 2 . 2 ) 


where the first determinant is taken over N v fields. The dependence on N is parametric, and the limit 
N —> 0 can be taken: 

, det (ilf — m v + J) „_s[A] 


ZReplica \J\ — 


\N V 


(2.3) 


det(iZ^ — m v 

It is not just that the two methods succeed in removing the quark determinant in the partition 
function. The two generating functionals ([O]) and (^]^) are simply identical. Passing to the effective 
field theory description, the equivalence is thus guaranteed from the very beginning. This can easily 


be verified at the perturbative level of the chiral Lagrangian with the same sources [11]. It is also 
clear that it trivially generalizes to the partially quenched case. 


We assume that in the quenched theories so constructed, the chiral flavour symmetry is sponta¬ 
neously broken to the diagonal subgroup. The low energy degrees of freedom are the corresponding 
Goldstone bosons, the dynamics of which can be described by a chiral Lagrangian. In the replica 
method, the symmetry breaking pattern is as in QCD SU(N)l x SU(N)r —► SU(N)r + r, while 
in the supersymmetric case we have naively a graded flavour symmetry of the breaking pattern 
SU(N v \N v )r x SU(N v \N v )r —> SU(N v \N v )l + r. One new feature common to both methods comes 
from the fact that in contrast with the full theory the flavour singlet cannot be integrated out: it is a 
degree of freedom that does not decouple. This is well-known in the supersymmetric formulation, and 
it is even more easily understood in the replica formalism. There is simply no replica limit N —> 0 of 
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a theory with SU(N) symmetry, and the trace part must be allowed to fluctuate]^ 


Although the perturbative series obtained by the two methods coincide, we need to go beyond per¬ 
turbation theory here, as we shall be treating the zero momentum modes of the Goldstone bosons 
exactly. Then a very precise definition of the supersymmetric Haar measure in the supersymmetric 
formulation, different from the one employed in formal perturbative expansions, is required. This has 
been explained in detail in refs. [ 12 , 


(see also the recent discussion 


To obtain the precise 


Haar measure over which to integrate, a more careful analysis of the flavour symmetries with bosonic 
and fermionic quark species is required. The outcome is that one should replace the naive supersym¬ 
metric generalization U(N v \N v ) by Gl(N v \N v ), or rather what in the mathematics literature is called 
the maximally symmetric Riemannian submanifold [12] of this supergroup (for simplicity of notation 
we will not make this distinction in what follows, and will just denote it by Gl(N v \N v )). It basically 
amounts to choosing well-defined integration paths for the fields involved; in particular assuring that 
the action is bounded from below in all bosonic directions. The replica method is in principle much 
easier since no supergroup is needed at all, and one is working with an ordinary chiral Lagrangian 
throughout. The perturbative Feynman rules are simpler, but at the non-perturbative level it is not 
known how to go beyond series expansions in general [ 151 (an exception is QCD 3 , where the expansion 
terminates |l(| and it thus gives the exact result [[17]]). So while equivalence is trivial in perturbation 
theory, it is highly non-trivial at the non-perturbative level. By doing all computations both ways, we 
thus get important cross-checks on our results. 


3 Regimes of (Quenched) Chiral Perturbation Theory 

Let us consider QCD in a toroidal volume V of average length scale L = R 1 / 4 . We assume that 
the volume is large with respect to the QCD scale. As in an infinite volume the lightest degrees of 
freedom are the Goldstone bosons of chiral symmetry breaking. They are describable in terms of a 
chiral Lagrangian given by an expansion in powers of the pion momentum p = ‘I'nn/L and mass m n 
over the cutoff of the effective theory, A ~ 4ttF. This is the standard chiral expansion, in which both 
quantities p and m n are taken to be of the same order. At leading order, the chiral Lagrangian in a 6 
vaccuum is usually writen as 

£ (2) = ^ Tr ^ dixUd[l U^ - ^TV (Me id/N fU + e~ id/N fU^M ] ) + O ^ , (3.4) 

where £ and F are the infinite volume quark condensate and pion decay constant, both in the chiral 
limit. j\4 is the quark mass matrix that for simplicity we will take as proportional to the identity, 
A4 = mF Furthermore, U = e l '^ 2 ^l F is an element of SU(Nf). In the chiral expansion U is expanded 
around the classical solution U = 1 + iy/2^/F + ... Consider the quadratic term in the action in 
momentum space 

= j(p 2 +ml) Tr[$ . (3.5) 

It provides a gaussian damping factor that limits the fluctuations of Tr[^] to be 

^ (yV(p 2 + ml)) O 

4 There is no corresponding problem with partial quenching, i.e. taking the replica limit IV —> 0 of a partially quenched 
theory of flavour group SU(Nf + N), with Nf > 2. Then the mo-term is not required to define the theory, and it can 
be decoupled. 
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When the linear size L is much larger than the Compton wavelength of the pions l/m T , the system 
hardly feels the finite volume, and the typical momentum scale is thus p ~ m n . So for m n 1/L, 
the fluctuations are much smaller than 1 / L 2 and the expansion around the classical solution is an 
expansion in powers of a quantity much smaller than 1 /(LF) 2 <C 1. As the size of the box or the 
quark mass becomes smaller, finite size effects start to become important, but provided m n L > 1 , 
ordinary perturbation theory is still applicable. In the boundary of this regime, when m^L ~ 1, the 
expansion of the field U(x) around the classical solution and the expansion in powers of momenta of 
the Lagrangian itself become the same expansion in powers of 1/FL. This is the so-called p-expansion 
|lj in which 

\ip\ p _. 1 

F ~ A ~ A ~ LF ' 


If the chiral limit is approached further in such a way that the Compton wavelength of the pion is 
larger than the box size (L > l/m n ), the conventional p-expansion eventually breaks down due to 
propagation of pions with zero momenta jl|]. Indeed, according to eq. (|3.5|) , when m^L 2 ~ 1 (in 
units of the cut-off, i.e. m 7T L 2 F = (2mSk) 1 / 2 ~ 1) the fluctuations of the zero modes, £q, are of 
0(1) and the perturbative expansion for these modes breaks down. But the non-zero modes still have 
fluctuations that are 



and they are still perturbative. A convenient expansion for this regime is the so-called e-expansion, in 
which [FT! 


my 

A 


v_ 

A 2 


1 


L 2 F 2 


The zero modes of the pion can be isolated by factorizing U(x) into a constant collective field JJq and 
the pion fluctuations £(x): 


U(x) 


U 0 exp i 


^(x) 

F 


(3.6) 


The difficulty in this regime comes from the fact that the integral over Uq needs to be done exactly, 
while ordinary chiral perturbation theory applies to the non-zero mode integration, since £ 2 rs_/ 0 (e 2 ). 
To leading order, the partition function is [] 


Z(6,M) = [ d^dUoexp 
Jsu(N f ) 


d A x Tr (3^3^) + ( e Wl7o + 


If FL 1, the p- and e-expansions should match in the range of quark masses such that m n L ~ 1 . 
In this regime rnFV ~ (FL) 2 (m n L) 2 ^$> 1 , so the results in the e-expansion should reproduce those of 
the p-expansion in the limit of large mT,V. 


It is also interesting to consider averages in sectors of fixed topology [|J. Fourier-transforming in 6, 
we get, to the same order: 


Z V (M) = / d£ [ 

JSU(Nf) J v 


d£ I dUo ( detUo) l/ exp 
SU(N f ) Ju(N f ) 


i J d A x Tr (d^O + —t(—T r (u 0 + t/ 0 f ) 


The e- and p-expansions can be defined in the same way in this truncation of the theory. 

s The measure contains additional terms of order e 2 and higher |:j[] due to the change of variables to coordinates £(*)• 
These terms are not needed to the order at which we will be working. 
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3.1 Quenched Chiral Perturbation Theory 


In the supersymmetric formulation of the quenched theory, assuming the pattern of chiral symmetry 
breaking of the previous section, the low-energy behaviour of the theory can be described by the 
supersymmetric chiral Lagrangian, 

£( 2 ) = ^Str (d^U~%u) - ^Sti (UgU + U- l U 0 l ) + ||<f> l + , (3-7) 

where Str denotes the supertrace, <f > o = -^=Str[— i log(C/)] and Ug = exp(i0 /N v )Ijy v + In v - Here, 
In v is the identity matrix in the fermion-fermion block of “physical” Goldstone bosons and zero 
otherwise, while Ijy v is the identity in the boson-boson block and zero elsewhere. As explained in 
the previous section, the integral over the Goldstone fields is over a submanifold of Gl(N v \N v ). An 
important difference with the full theory is the unavoidable presence of the singlet field, which cannot 
be decoupled in this case. Actually, the chiral expansion signals the breakdown of the perturbative 
series if the singlet mass mo becomes comparable to or larger than the cutoff scale in the effective 
theory A ~ 4 nF. To see this, let us consider the case of N v = 1 in the context of the supersymmetric 
method. Using the parametrization of the Goldstone submanifold of G7(l|l) from ref. |13|], 


U(x) 



0 


7 (*)^ 

0 


where 4>(x) and s(.x) are the real bosonic fields, and ^(x) and (3(x) are the fermionic ones, we consider 
the usual p-expansion about the classical vacuum U = 1. The action expanded to quadratic order is 



where m^F 2 = 2mS. In contrast to the naive choice of 17(111) as the coset of chiral symmetry breaking, 
the gaussian integrals are absolutely convergent for all momenta such that (mg + ap 2 )/N c < p 2 + m 2 
and to any finite order in (ttiq + ap 2 )/N c even if this condition is not satisfied. It is only then that we 
are able to make meaningful statements about the magnitude of fluctuations. 

Using the properly normalized measure we get to this order: 

_ det(-<9 2 + ml) 

det(— d 2 + ml) 1 / 2 det(— d 2 + ml) 1 / 2 

and Z is thus independent of mg in the absence of sources. This explains, to this order, the puzzle of 
how the effective partition function apparently could depend on one new parameter mo, when, in the 
absence of sources, it should equal unity by construction. Only when we include external sources for 
the quarks can a dependence on mo appear; in particular individual propagators will depend on mo¬ 
lt follows that mo must have at its origin gluonic dynamics, which is probed by appropriate quark 
sources. This picture is in nice agreement with what is found at the level of the zero modes, see below. 


(3.9) 


The fluctuations of the Fourier modes of s and 4> are seen to be 



(1 - a/N c )p 2 + ml - ml/N c 
V(p 2 + ml) 2 



(1 + a/N c )p 2 + ml + ml/N c 
V(p 2 + ml) 2 


(3.10) 
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Clearly if rriQ/N c is larger than the chiral scales p 2 or m 2 , the fluctuations of the fields are controlled 
by this parameter and the perturbative expansion breaks down for all modes (i.e. zero and non-zero) 
when mg/iV c > F 2 []. This is an artefact of the quenched approximation, related to the fact that the 
singlet cannot be integrated out. Defining e /2 = (m^)/ (N C F 2 ) and e 2 = 1/(FL ) 2 , in the regime of the 
p-expansion we have 


1 

WW 


9 

rn * 
A 2 


0(e 2 




~ 0(e 2 + e /2 ) . 


(3.11) 


In the case of the zero modes, eqs. ( |3.10| ) imply that the series breaks down when m 2 / A 2 < e'e 2 . For 
e' > e 2 this will happen before zero modes need to be treated separately in the unquenched theory. 
This is, however, not the case if we only consider averages in sectors of fixed topology at finite v. Using 
left and right invariance of the super-Haar measure, the factorization into zero modes and non-zero 
modes goes through as in the unquenched case eq. ( |3.6| ), and the partition function restricted to fixed 
topology v is to the order we need it 


MM) = 


o-v 2 !^ 2 ) 


yj2ir(i 


/ Gl(Nf\Nf) 


dUod £ (SdetI/o)^exp 


'mTiV 


Str(t/ 0 + M 1 ) 


+ jd 4 x(~Sti(d^dM - ^-Str(£) 2 - ^ (0 M Str[£]) 2 ) + ... 


where Sdet is the superdeterminant. The distribution of winding numbers is Gaussian with mean 

F 2 mQV 


(O = 


2 Nr . 


(3.12) 


a relation that will play a crucial role when we match results in the e-regime with those of the p- 
expansion regime. The otherwise menacing <ho has disappeared in the zero-mode sector at fixed 
topology, leaving as its only trace the average distribution of topological charges. As a result of this, 
the usual e-regime, where 


1 


(FLY 


m n ✓-»/ 2 ^ 

~ ~r ~ °( e ) 


A 


(3.13) 


can be approached. This result may seem puzzling, as it indicates that a more chiral regime can be 
reached in sectors of fixed gauge field topology, whereas after summing over topology we find that 
the double-pole term of the quenched propagator blows up at a scale much before the regime of the 
e-expansion is reached. In fact this is correct, and the resolution is found by noting that also the 
quenched e-expansion in sectors of fixed topology eventually fails, when \v\ —> oo. In the sum over 
topology, the dominant contributions are those around \ v\ ~ Fitiq^/V and for topological charges that 
large, the e-expansion breaks down even in sectors of fixed topology. In other words, the sickness that 
we found in the direct analysis for m 2 /A 2 < e'e 2 reappears in this case when summing over topology. 
Actually this problem of the perturbative expansion close to the chiral limit is probably unresolvable, 
since it is also found at the quark level: in a background gauge field with topological charge u, the 


contribution of the topological zero modes to any quark propagator is \v\/mV. Using eq. ( 3.12 ) the 
average over topology implies that this contribution becomes of order Fmo/{my/V), which is of 0(1) 
when m ~ e'e 2 . Concerning the non-zero modes, the perturbative expansion remains a simultaneous 
expansion in e and e' also at fixed topology. We will be working to lowest non-trivial order in both 
expansion parameters. 


“We assume the parameter a ~ 1. 
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Let us now summarize the relevant Feynman rules. We have already displayed the usual supersym¬ 
metric version of the Lagrangian and the propagators for the bosonic fields in eq. ( |3.1C ) for the N v = 1 
case. The propagator for the fermionic fields 7 and /3 is, after a simple rescaling, read off from eq. 
( |3.8| ) to be simply a conventional bosonic propagator. These rules differ by some signs from the usual 
ones based on 1/(1|1) |l(J. These differences are irrelevant in perturbation theory. The generalization 
to N v > 1 follows similarly [|l(]], with a few sign changes. Let us introduce some convenient notation 
and define 

1 p ipx 1 p ipx 

A M = fEro . A M = y'E' — • < 3 - 14 ) 


V ^ p 2 + m 2 


P 


where a prime on the sum indicates that zero momentum is excluded. We will also need the corre¬ 
sponding expressions for the peculiar double-pole term. Let 


G{x) 


1 e v x 

- y p {p 2 + ml) 2 


and G(x) 


1 p ipx 

_ _ / __ 

= 7 ^ p 4 ' 


One notes that 


d 4 x A(x) = / d^x G(x) = 0 


(3.15) 


(3.16) 


properties that greatly facilitate comparison with integrated Ward identities in the e-expansion. There 
is also the relation 


A(x) = 


1 


mzV 


+ A(x) - mlG(x) + ... , 


(3.17) 


so that both A(x) and G(x) immediately follow from Taylor expanding the finite-volume massive pion 
propagator A(x). 


In the replica formulation [|ll]], all of the above discussion also applies. The action is that of eq. (3.7) 
after replacing Str by Tr everywhere, and after changing the group manifold to U(N). The mo term 
actually serves to interpolate between SU(N ) and U(N): when mo —> 00 we go from U(N) to SU(N), 
but then there can be no replica limit. The way the mo term serves to allow a replica limit becomes 
particularly transparent when one looks at the Feynman rules: in a quark basis the propagator of the 
off-diagonal mesons will have the usual form of A(x), while the diagonal combination has a propagator 
(with E being a N x N matrix with unity in every entry): 


1 


Gij ( X ) — 


J i3 


(V 2 + ml) 


- E, 


C ml + ap 2 )/N c 

lJ (p 2 + ml) 2 E(p 2 ) 


a ipx 


Here 


Hp z ) = 1 + 


mg + ap~ 

K 


N 


p 2 + ml 


(3.18) 


(3.19) 


For any finite N we can take the limit mo —■► 00 , in which case the double pole in the last term is 
cancelled. We then get the ordinary propagator with a factor of 1/N in front, combining with the first 
term to give the usual diagonal propagator in the quark basis. The singularity at N = 0 is regularized 
by keeping a finite mo- The replica limit N —> 0 is then taken only at the end of the calculation, 
keeping mo finite. 

4 The Chiral Condensate 


The leading order result for the chiral condensate in a finite volume and fixed topological sector was 
computed from the fully and partially quenched chiral Lagrangian in ref. [1^]- For full QCD the first 
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correction to this result follows directly from the calculation of Gasser and Leutwyler [1J. It is also of 
interest to find this first correction in the fully quenched theory, and this was recently done using the 
replica method [§]. In this section we show how the same result is obtained from the supersymmetric 
method. 


Let us consider the simplest case of N v = 1, because the result clearly does not depend on N v . For 
the group G7(l|l), the chiral condensate in a sector with fixed topology can be defined as: 


^(™) = yjJ h ' Z v( Mj )\j=0 ' 

where M.j = diag(m + j, m). At leading order in the e-expansion one finds [^, [i~8f 

M {IMKM + I„ + 1 {ii)K u _ i(/i)) + " 

where I u (n), K u (fj,) are modified Bessel functions, and /i = mEV. 




(4.20) 


(4.21) 


To derive the first correction to this result in the e-expansion it is convenient to first calculate the 
1-loop improvement of the action due to the fluctuations of the non-zero momentum modes. To 0(e 2 ) 
the contribution from the measure does not affect the result, and to that order we simply evaluate 

E 


1 - 


2 F 2 


Str 


M.J (Vo + t/ ( 7‘) J d A xe{x) 


Performing the integral over £ to leading order, and then re-exponentiating the correction, we obtain 

[ZeffV. 

'jKJ (J ^UUCUU(J ) — 

/Gi(l|l) 

Then the effective coefficient determining the strength of the condensate at finite volume is 


MMj) = J dU 0 (SdetC/orexp ^^Str (Mj (u 0 + U^ 1 )) 


(4.22) 


S eff(V) = S 1 + 


1 


Nr.F 2 


(mgG(O) + oA(0)) 


(4.23) 


were G, A are defined in eqs. ( 3.14 - 3.15 ). The mg and a terms are what remains from the partial 
cancellation of the fermionic propagator and the bosonic ones. The partition function is then the same 
as the one at leading order with the change E —> E e ff(V). 

It is now easy to calculate the condensate to one loop, by differentiating eq. ( 4.22 ) as in eq. ( 4.2C|) to 
get 

Si- lo °P( / r) = E v (//') — 




1 


= E^(/i) + 2^-^(m5G(0) + aA(0))g4(/i)A;(g) + 


where 


/r' = mYi e ffV = fj, I 1 + 


NrF 2 


(mgG(O) + aA(0)) 


(4.24) 

(4.25) 


One peculiarity of the quenched approximation is that the 0(e 2 ) correction to the condensate is 
ultraviolet divergent even in dimensional regularization because of the double pole propagator Jl^, |2^] . 
We find || 

G(0) = ft + ijHL/LoJ-d) , A(0) = -§, 
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(4.26) 



















where f3\ and /?2 are two of the universal “shape coefficients”. It depends on the shape of the box, 
and the precise value can, for any given volume, be computed from the general expression given in 
ref. [:5|]. In dimensional regularization the constant c\ reads 


ci = ^^-i(r / (l) + l + ln(47r))+0(d-4) 


(4.27) 


and 1 /Lq is the ultraviolet subtraction point. The divergent term in ci matches precisely the 1- 


loop counterterm needed for the condensate in the infinite-volume theory [20]. Below we present the 


numerical results corresponding to the subtractions that are required by the MS scheme. 


Note that the one-loop contribution contains a quenched finite-volume logarithm [j!|], which reflects 
the infrared sickness of the quenched approximation. It is the finite-volume counterpart of the usual 
logarithmic chiral divergence in infinite-volume quenched chiral perturbation theory. The analysis is 
clearly restricted to domains where this divergence does not yet overwhelm the tree-level result. 


Several methods have been proposed and implemented to calculate the coupling E by finite-size scaling 


techniques. One method is to take the discontinuity of the chiral condensate (4.24) across the imaginary 
jU-axis to get the spectral density of the Dirac operator eigenvalues in this regime [13]. Since the 
density depends on the parameter E, the leading tree-level expression for the density can be used to 
extract E in this way |n]]. Alternatively, the condensate may also be computed directly on the lattice 
|2^ . 24, |^] and one can then fit its finite size and quark-mass dependence to the prediction in 

eq. ( ^1.24|). B oth methods have been shown to work well in extracting E by neglecting the contributions 
of eq. ( |4.24|) . With sufficiently high statistics, one could also aim at a determination of the pion decay 
constant F [jl], |j] from these corrections. In the quenched case the problem is that it brings in two 
new unknowns: mo and a. To estimate the typical size of the correction, we can compare the 1- 
loop-improved chiral condensate at two different volumes, after subtracting the trivial topological zero 
mode contributions and keeping /r = mEk fixed. In this way we do not have to address the issue of 
the Lq scale dependence^]: 


z'drM - mm 

£;;i" p (A<) - M/m 


= 1 + 


1 1 LiI u (n)K u (n) 

4tt 2 N c F 2 S^/x) - M lii 


m » M (h) _ oft (4j - T 


(4.28) 


As an example, for the smallest and largest lattice volumes considered in [23] in the v = 1 sector, the 
second term in eq. ( 4.28 ) is around 15% for mo ~ 600 MeV. Although the statistical errors in [26] are 
larger than this, this is not clearly beyond hope in more precise studies. 


5 Quenched Correlation Functions 

Meson correlation functions can be computed systematically in the e-expansion. Here we present 
the results of the 0(e 2 ) calculation for scalar and pseudoscalar correlation functions in the quenched 
theory. We have performed all computations in both the supersymmetric and replica formulations, 
checking that the final results are identical. We begin with the flavour-singlet combinations involving 


S (x) = tp(x)I Nv ip(x) , P u (x) = 'tp{x)i^ b I n v iP(x) , 


(5.29) 


where 7 jv„ is the flavour projector onto the physical quark sector or that of the N v valence quarks. For 
simplicity we consider N v = 1 for the singlet case. To compute these correlation functions, we make 

7 In the direct comput ation of the quark condensate UV divergences appear at the quark level which must be subtracted 
before comparing to eq. (4.28) ||2:j(]. 
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use of the correspondence between local sources added at the quark level 

£qcd —> Cqcd + s(x)S°(x) + p(x)P°(x) , (5.30) 

and the substitution 

M ->• x(x) = M + s(x)I Nv + ip(x)I Nv (5.31) 

in the effective theory. Two functional derivatives of the generating functional with respect to s(x) 
and p(x) thus give us the scalar and the pseudoscalar 2-point functions, respectively, to the desired 
order. In the supersymmetric formulation we follow the procedure outlined in section 2, and thus add 
the sources in the fermion-fermion block only. Using the two sets of Feynman rules for the fluctuation 
fields £(x), we find, to order e 2 : 


(S°(x)S»> 

(P°(x)P°(0)) 


with 


c°p- 


X 2 
2 F 2 
X 2 
2 F 2 


(™%G(x) + qA(i)) a_ - A(x) a+ +a 2 -- 

( m oG{x) + «A(s)) «+ - A(x) a+ + ^~ +4 


(5.32) 

(5.33) 


a ± = <(^ii ± ^i 1 ) 2 ) • 

The constant terms are given by the same expectation values of eq.(5.34) 


X 2 

riO _ Zj e//„l~^oop 


Cl = - 


X 2 

2 ~‘ e ff „1 -loop 

4 °- 


(5.34) 


(5.35) 


but now evaluated with respect to the one-loop-improved action ( 4.22 ), and X e jy is as defined in eq. 
(4.23). In this way we obtain these constant terms to the required order e 2 . 


These expressions are valid in both formulations. As a first check, we note that when U —> 1 the 
x-dependent part of (S' o (x)S’°(0)) vanishes: there is no tree-level scalar propagation in the quenched 
theory. But there should be tree-level propagation of the quenched pseudoscalar singlet, and indeed 
in this limit the x-dependent part of (P°(x)P°(0)) approaches 


2X 2 



t G[x) 


(5.36) 


which is simply the massless singlet propagator (cf. eq. ( 3.18 )). 
comparison, also including the zero-momentum modes, below. 


We shall make a more detailed 


To evaluate the remaining group integrals in closed form, we can make use of the exact results obtained 
in the supersymmetric formulation |jl3|. For a source /xj = p + J in the fermion-fermion slot, the 
generating function becomes 

Z V [J] = ^pj (I u+ i(pj) + I u -i(pj)) K v {p) + ^pl u (pj) (K u+1 (p) + K v -i{p)) . (5.37) 

Using the fact that Z v [ 0] = 1 (which here is a consequence of two Wronskian identities for the Bessel 
functions), this gives us the expectation value 

((Un + (t/- 1 )!!) 2 } = 4 £pZ v [J\ 
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= 4 


— Iv+1 i(aO + 1 + 


K^-i) 


/+ 


(5.38) 


It is quite simple to find the remaining matrix element by means of the explicit parametrization: 

( UuiU -^ u ) = (1 + 7/3) 

r2n slQ roo 

= 1 + 2 / — dsd'yd/3 1 f 3e ^e-s) e ^cos{e)-cosh(s)) 

Jo 27t Jo 

= 1 + 2 I u (p)K u (p) . (5.39) 


These expectation values suffice to compute the correlation functions in eqs. ( 5.32 ) and ( 5.33| ). Rewrit¬ 
ing the final answer in terms of the tree-level chiral condensate eq. (4.21) and its derivative 


^ = I u (p)K u (p)-I u+1 (p)K u _ i(m)-4> 

Zj fl 


we Hnd: 


a+ = 4 


Kfa) v 


+ i + 




, a_ = 4 


1 E„(/i) + u 2 


L 




(5.40) 

(5.41) 


while the constants to the same order are: 


r° - 

l '5 — 


= 


J e// 


V '2 




+ 1 + 


J e// 


1 £„(//) ZZ 


= 5+ 


s^ oop '(M) 


+ 


¥)' 


+ 




,/2 


= E 2 

'l T}~ lo °P{p) v 2 ' 


p y p? 


y ^eff F' 

where p! and E e py are defined in eqs. ( 4.25 ) and ( 4.23|) , respectively. 
Using the relation 

y/ („:\ ^ y,\ 

as p 0 , 


(5.42) 


K(l) i s y (/x) 

S /r S 


0 


(5.43) 


it is easy to show that in the sum of the scalar and pseudoscalar correlation functions, the poles in 
the quark mass due to the zero modes cancel, and the sum has then a well-defined massless limit even 
at finite volume for v ^ 0. The case v = 0 is special: there is then an additional infrared singularity 
due to the distribution of the smallest non-zero Dirac operator eigenvalue. However, the infrared 
divergence for p —* 0 in that case is only logarithmic. 

It is instructive to see, analytically, how the p- and e-expansions match one another when m n L ~ 1 
(with FL 1). We have already noted above that the x-dependent part of the (P°(x)P°(0)) correla¬ 
tion function to 0(e 2 ) approaches the tree-level propagator of a massless flavour singlet. Now that we 
have the constant part of this correlation function evaluated explicitly, we can see how this constant 
part precisely serves, to leading order, to restore the full massive propagator of the p-expansion. To see 
the matching we need to sum over topology first, because this is how the p-expansion is usually done. 
After performing this sum over topology, using relation (3.12) as well as the Gell-Mann-Oakes-Renner 
relation, we get, to leading order: 


(P°(x)P°(0)) 



1 - 


a 


A(x) - ^G(r) 
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2 E 2 

~W 

2 £ 2 

“F2" 


a(i) -v;£ 


(y->2 pipX YYl 2 

PT^p-^ G(l ) + 0(m ’) 


< 7 (x) + 0 (m 2 ) . 


(5.44) 


As in the case of the one-loop correction to the chiral condensate the sum over topology precisely 
restores the mo-dependence needed for the e- and p-expansions to match. 


For comparisons with lattice gauge theory data, it is convenient to recast the above results in terms 
of the space-integrated correlators in the (euclidean) time direction, here labelled by t: 


V\t) = / d 3 x{P\x)P\ 0)) 


S\t) = / d 3 x{S\x)S\ 0)). 


(5.45) 


For this purpose we need the projections onto zero momentum of both the usual massless propagator 
and the double-pole term. They both follow immediately from Taylor-expanding the massive pion 
propagator in powers of m n . From 


Jd 3 x A(x) = 


cosh(m ff (T/2 — t )) 
2 m 7 r sinh(m 7 r T/2) 


and eq. (|3.17l) , we find 

/Aa(x) = ^ + |[( T -i) 


1 

12 


+ 


T 3 

24 


r 2 (r-l) 2 - 


30 


ml + 


which makes it convenient to define 

hi{r) = ^ 


d 3 x A(x) = 


1 


M r ) = -^3 J d 3 x G{x) = 


(r- -) 2 - — 

v 2 J 12 


1 


24 


r 2 (r-l) 2 - 


30 


(5.46) 

(5.47) 

(5.48) 

(5.49) 


Here T is the extent in the temporal direction, and r = t/T. These 2nd and 4th order polynomials 
replace the exponentially decaying correlations in the regime of the e-expansion. 


Figures |l] show the result for P°(t ) and S°(t ) in this e regime, for three values of the topological charge 
v = 0,1, 2, and a quark mass m = 3 MeV in a typical lattice of V = 16 3 32 at a~ 1 = 2 GeV. In these 
plots, and in those of the next sections, we have used the values mo = 600 MeV and a = 0.6. For £ we 
take as indicative the value (270 MeV ) 3 at the commonly used renormalization scale Lq = ( 4:irF ) -1 . 


We should stress that in these predictions, and in any other predictions given in this paper, the contri¬ 
butions from higher excited states with the same quantum numbers are neglected. These higher states 
clearly do contribute to the correlation functions at short time separations, but they are exponentially 
suppressed at large t. Therefore our predictions should only fit the measured correlation functions for 
sufficiently large time separation between the sources. 


5.1 Flavoured Correlation Functions 

The computation of mesonic correlation functions in non-singlet channels is a little more difficult, 
because we need the exact evaluation of the group integral over zero momentum modes for N v > 1. 
Here we restrict ourselves to the simplest case of N v = 2. 
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Figure 1: pV°{t) (left) and pS°(t) (right) in a lattice of V = L 3 T = 16 3 32 and a 1 ~ 2 GeV. The 
quark mass is m = 3 MeV and the three curves correspond to v = 0,1 and 2. 


Let us denote the quark bilinears as follows: 

S a (x) = ip(x)t a I Nv ip(x) , P a (x) = 'ip(x)t a i'y 5 I Nv 'i/;(x ) , (5.50) 

where in this case t a = r a / 2, and the t“’ s are the usual Pauli matrices. The corresponding two- 
point functions are (S a (x)S a ( 0)) and (P a (x)P a ( 0)), for a=l-3. Here we present some details of the 
calculation of the correlators for a = 1. 


The appropriate flavoured sources s a (x ) and p a (x ) are obtained by lettting 

£qcd -> £qgd + s a (x)S a (x) + p a (x).P a (x) 


(5.51) 


at the quark level. This translates into 


M -»• %(x) = X + s a (x)t a I Nv + ip a (x)t a I Nv 


(5.52) 


in the effective theory. Two functional derivatives with respect to these sources then give us the 
desired correlation functions. To order e 2 , we find: 


(S\x)S\0)) 
(P\x)P\ 0)) 


C 5 + 

Cp- 


s 2 

2F 2 

S 2 

2F2 


^mgG(x) + aA(x c_ — A(x)b_ 
(m§G(x) + aA(i)) c+ - A(x)6+ , 


where, after using some of the identities provided in Appendix |A[ one finds 

c± = - ({U 12 + U 21 ± (17 X )l 2 ± (t^^ 1 ) 2 l) 2 ^ 
6± = ^(c/iiC/ 2 2 + (C/“ 1 )n(^" 1 )22)±l • 


(5.53) 

(5.54) 
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(5.55) 

(5.56) 
















The constant terms will also be evaluated to 0(t 2 ). Since they are flavour-symmetric, we can simply 
take the constant part of 

1 z , n 

3<S»“(i)<Ss“(0) 1,1 1 

and similarly for Cp. Using the SU( 2) completeness relation 


j a j .a 




(5.57) 


we then get 

C s = ((Un + (U- 1 ) n) 2 + (i U 22 + (U- 1 )^ 2 ) + 2(Un + (U- 1 ) 12 )(U 21 + (U~ l ) 2l ) 

- (Un + (U~ 1 )ii)(U 2 2 + (U -1 ) 22 )) 

C P = ((^11 - (^- 1 ) n ) 2 + (U 22 - (U~ l ) 22 ) 2 ) + 2 (U 12 - (U~ l ) 12 )(U 2 i - (U~ l ) 2l ) 

- (Un - (U _1 )n )(U 22 - (U-% 2 )) . (5.58) 


Here the expectation values are again to be computed with respect to the one-loop-improved action 
eq. (4.22) so as to reach 0(e 2 ) accuracy. 


In order to calculate the various integrals over the zero-momentum modes explicitly, we make use of 
the fact that they can all be related to an expectation value already computed in eq. ( 5.38| ), and 


((Un + (U- 1 )n)(U 22 + (U- l ) 22 )) 


(5.59) 


which was evaluated explicitly for the zero-momentum mode integral with N v = 2 in ref. |27|]. The 
idea of how to get the other integrals from these two is simple, and most easily explained in the replica 
formalism. There it is particularly clear that the quantity det(JUl)~ u Z v is a function of the eigenvalues 
of A4A4^ only. This means that all quadratic expectation values, after appropriate manipulations, can 
be related directly to the expectation value of eq.( 5.59|) , and that of eq. ( 5.38 ). 


Using a variety of 2 x 2 matrix sources J, we have collected a list with the results of explicit evaluations 
of all needed zero-momentum mode integrals. For the reader interested in the technical details, see 
Appendix [A[ Here we simply present the final results: 


and 


b+ 


c+ 


2 ^1 + 




C- 



2 1 ) 
fJL S 


= p 2 ^y oop V), 


Cp = . 

2 /x 


(5.60) 

(5.61) 


As a non-trivial check on the calculation, we have also evaluated the correlators corresponding to 
a = 2 and 3. In the latter case the intermediate details are very different, but the final results coincide 
in both cases with the ones above for a = 1 , as expected. 
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In the sum of the two flavoured correlation functions, our expressions simplify. We find (no sum on 
a): 


(S a (x)S a ( 0)) + (P a (x)P a (0)) 




A(x) — 



+ ^ 5 E^)j^( m ” eW + “ AW ) 


(5.62) 


As in the flavour-singlet case, the sum remains finite in the p —> 0 limit for v ^ 0, on account of eq. 

(ED- 


Here we can also check analytically that the e- and p-expansions join on to each other. This is because 
the (P a (x)P b ( 0)) correlation function we have computed should simply reproduce the tree-level pion 
propagation of the p-expansion to 0(e 2 ). Indeed, we find, to lowest order (no sum on a): 


(P a (x)P a { 0)) ~ 


£ 2 S 2 - . 

v + i^ A(x 

S 2 


F 2 

£ 2 


1 


miV 


+ A(x) 


A(*) + °( m l) 


(5.63) 


which is precisely the tree-level result of the p-expansion. The constant again provides the term needed 
to restore the pion zero-mode term in the propagator. 


Again it is convenient to reformulate these results in terms of projections onto zero momentum. The 
relevant definitions have already been given above, and in Figure || we show V a (t) for u = 0,1, 2 and 
various quark masses within the e regime. As expected, the effect of the zero modes are dominant in 
this regime, inducing a large dependence on the topological charge. The same effect is observed for 
the scalar correlator in Fig. [3]. Note that this correlation function is negative for v = 1,2. This feature 
was first noted by the authors of [31] for larger quark masses in the p-expansion regime (ie. not at 
fixed topology). As explained above, it is possible to cancel the poles induced by the topological zero 
modes by considering the sum of the scalar and pseudoscalar propagators. This is shown in Figures ||. 
The only divergence in the p —> 0 limit is then restricted to the v = 0 sector and it is only logarithmic. 


As far as we are aware, there are at present no lattice data for any of these correlation functions in 
quenched QCD at fixed topology and in the appropriate quark mass range. Very recently, there has 
been a study of pseudoscalar correlation functions in the Nf = 2 Schwinger model at fixed topology 
p^| , and distinct signatures of topology were actually found there. But, unfortunately, neither the 
present calculation for quenched QCD nor the one presented in section |6] for full QCD at fixed topology 
can directly be compared with these two-dimensional results. It has also recently been suggested to 
study QCD at fixed topology, but at volumes so large that correlations are still exponential in the 
masses |2^[ . 


5.2 Ward Identities at Fixed Topology 

Ward identities associated with the chiral rotation of the physical (valence in the replica context) 
quarks provide useful cross-checks on the above results. For the singlet chiral transformation, the 
relation is 


{(d fl A°(x)-2rnP°(x)-2iN v uj(x))O(0)) = - (5O(0))5(x), (5.64) 
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Figure 3: fiS a (t) for a quark mass m = 3 MeV and u = 0,1, 2. The lattice volume is V = L 3 T = 16 3 32 
and a -1 ~ 2 GeV. 




Figure 4: V a (t) + S a (t ) for quark masses m = 0.1,1 and 3 MeV and v = 0 and 1. The lattice volume 
is V = L 3 T = 16 3 32 and a -1 ~ 2 GeV. 






(5.65) 


where A ° = i/j!jv v 7 /j7sV’- The operator 


cj(x) = 


lg7r 2 TrF ^ F ^ 


is the topological charge density, 0(x) is any local operator and 50(x) is its variation under an 
infinitesimal singlet chiral rotation at point x. The identity of eq. (|5.64| ) is normally considered in 
the full theory, after summing over all topological sectors. But it is exact, and in fact simpler, also 
in sectors of fixed topology. We also remark that although all effects of quark loops disappear in 
the quenched theory, the above identity reflects the consequences of a chiral rotation in part of the 
theory only (i.e. in the valence sector). It would reduce to a triviality if we were to rotate all fields 
simultaneously. 

If we combine the Ward identity for 0(x) = P°(x) and 0(x) = uj(x), and integrate over space-time, 
we can eliminate (uj(x)P°( 0)) to get at fixed topology 

/ 2 (S°) 


Jd 4 x (P°(x)P°( 0)) = - 


m 2 V mV 


(5.66) 


where S° is defined in eq. ( 5.29|) and we have here used N v = 1. Inserting eq. (5.33) and taking into 
account eq. (3.16), we get, to 0(e 2 ): 


C° P = 


K~ loop {n) v 2 


mV 


(mV) 2 


(5.67) 


which coincides precisely with our result in eq. ( |5.42| ). Actually, at fixed gauge field topology, 
(cn(x’)P°(0)) can be computed directly using eq. QA.17| ) of Appendix |Aj: 

i<p») = VL . (5.68) 

/X 


Inserting this into the identity of O(x) = P°(x) gives precisely the same relation as ( |5.67p . 

Similarly, the general Ward identity for a flavour non-singlet chiral rotation is 

{(d, l A“(x)-2mP a (x))O(0)) = -(5 a O(0))5(x) , 


(5.69) 


where 5 a O(x) is the variation of O under an infitesimal chiral rotation at x in the flavour direction 
a. Combining the identities thus obtained for 0(x ) = P a (x ) and 0(x ) = <9 M A“(x), and subsequently 
integrating over space-time we arrive at the relation, to order e 2 (no sum on a): 


Jd A x (P a (x)P a ( 0)) = - 


1 (5°) 1 Y} u - loop {v) 


4 mV 


+ 


m 


(5.70) 


which also agrees with the constant term Cp of eq. (5.61), computed to the same order. Note that S° 


in eq. ( 5.70 ) is the one corresponding to N v = 2. 


Although not Ward identities as such, a series of relations that follow from spectral representations 
of the involved correlation functions can also give useful checks on our results. Such relations were 
derived in ref. [f30| , using the language of lattice overlap fermions, but they can easily be transcribed 
into continuum language. In the quenched case a particularly interesting relation gives a sum rule for 
the flavoured scalar correlation function. In our normalization, with no sum on a, 


r d A x (S a (x)S a ( 0)) = - 477^-(S 0 } = 

\ \ J \ J/ 2 V dm x ' ° 


+ 


(5.71) 
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which agrees with our value of the constant in eq. ( 5.61 ). The authors of ref. [30] have also 
derived some general relations for the m —> 0 limit of the finite volume correlation functions we have 
considered here. For N v = 2 they are in our normalization (again no sum on a): 



(P°(x)P°(0)) — 4(S a (x) S a (0)) 




Jd A x [(P a (x)P a ( 0)) - (S a (x)S a (0))] ~ 


4|i/| 

m 2 V 

|i/| 

m 2 V 


4 v 2 
m 2 V 


(5.72) 

(5.73) 


as m —> 0. We have not computed (P°(x)P°( 0)) explicitly for N v = 2, but to check eq. ( |5.72|) we 
need only the constant part of this N v = 2 correlator. It is fixed by the Ward identity of eq. (|5.64[) to 
be 

N N v u 2 

v mV (mV) 2 


to lowest order. Plugging this in, and using the relation eq. ( 5.71 ), which we have just confirmed, we 
find that eq. ( 5.72| ) is indeed satisfied. Similarly, eq. ( 5.75 ) is easily seen to be satisfied by taking the 
m —> 0 limit of our exact finite-volume expressions. 


6 Full QCD at Fixed Topology 


Eventually lattice computations with light quarks will be pushed beyond the quenched limit. It is 
therefore useful to know the corresponding analytical expressions for finite-volume meson propagators 
also in full QCD. The e-expansion of the chiral Lagrangian is much simpler in this case, as all effects 
of the flavour singlet term can be ignored. More precisely, we obtain the relevant chiral Lagrangian 
by taking the limit mo —> oo of the replica chiral Lagrangian. This enforces Tr£(x) = 0, and the 
singlet field decouples. The relevant expansion was carried through to order e 4 in ref. Q, probably 
beyond the realistic accuracy of lattice simulations in the near future. It is of interest to find also the 
corresponding analytical expressions in sectors of fixed gauge field topology. Such expressions cannot 
immediately be inferred from ref. Q, as all results are listed in a manner valid only for flavour group 
SU(Nf) 8 9 . We first list the correlation functions up to order e 2 in the full theory with summation over 
topological charges, which agree with the results of ref. || to this order[]: 


<S°(z)S») 

(P°(*)P») 

(S“(i)S»(0)) 

(P“(z)P‘{0)) 


y 2 

r° — 

C s 2F 2 


C°p + 


P ' 2F 2 

5 ab r 


(Tr [(U - Dt) 2 ]) - 1- ((TV (U - t/t)) 2 ) 
(Tr[(t/ + Dt) 2 ])-^((Tr(C/ + t/t)) 


A(a;) 

2 ^ A(x) 


N f~ 1 


C S - 


4 F 2 


((Tr(D)) 2 + (TV([7t)) 5 


Nj 

5 ab 


2 

(Tr (U 2 ) + TV(D t2 )) - 2Nj + 4 - ^ ( Tr(U)Tr(U f ) 


N f~ 1 


C P + 


4 F 2 


AM) ((TV(l7)f + fIY(t/t)) 2 


Nj 


A(x) 


8 In particular, the author of ref. Q has implicitly used the identity (Tr(C/)) = (Tr(Lb)}, which holds only after 
summing over topological charges. Also, other identities specific to the group SU(Nf) have been employed prior to the 
listing of results in ref. Q. 

9 Our notation is the same as in the previous sections. 
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- (Tr (U 2 ) + Tr(tf 2 )) + 2N] - 4 + A (Tr([/)Tr(t/t)} 

f f 


A(x) } . (6.74) 


To the order at which we are working, we need the 0(e 2 ) contributions to the constant terms. As 
in the quenched theory, these are entirely given by the one-loop correction to the chiral condensate, 
which here reads [ill 


Xeff (V) -I , N f ~~ 1 1 Pi (Li/L) 
£ N f F 2 L 2 


(6.75) 


Here j3\(Li/L) is another of the universal shape coefficients J2], [5| . To order e 2 the constant terms are 
given by 


eg 

C% 


C s 


Cp 


^"((Tr([/ + C/t)) 2 } 
-^"((T7(C/-C/t)) 2 } 
^ (Tr[(C/ + C/t) 2 ]}- 



(Tr[(t/-C/t) 2 ]) 




(6.76) 


where the expectation values are taken with respect to the one-loop improved action (simply let 
E —► in the tree level action). As a check on these results, we note that upon substituting U —* 1 

the coefficients of A(x) vanish for the first three correlators. Indeed, there should be no tree level 
propagation of non-zero modes in any of these correlation functions 10 . A similar substitution in the 
last (flavoured, pseudoscalar) correlation function yields the coefficient E 2 /F 2 , corresponding to tree 
level propagation of the non-zero-momentum modes in this channel. 


Next, to project down on sectors of fixed topological charge v we again integrate over 0 as in section 
H|, simultaneously extending the group integration from SU(Nf) to U(Nf) for the zero-momentum 
modes. Thus only the integration over the zero-momentum modes is affected by this projection, and 
the fluctuation part is precisely as in ref. ||. In particular, the fluctuations £ a (x) still belong to 
the adjoint representation of SU(Nf). The above expressions then remain valid in sectors of fixed 
topological charge, once expectation values over the zero-momentum modes are interpreted as being 
with respect to the measure of 

Z v = f dU (det U) u exp 
Ju(N f ) 

One obvious benefit of going to sectors of fixed topological charge is that all pertinent expectation 
values can be evaluated explicitly for any Nf. In order to do so, we first invoke the identities (see 
Appendix A for a derivation) 


^Tr(U + Uy 


(6.77) 


(Tr([/ 2 )) = N f - 2 ^ Nf + ^ (Tr(17)) 

(Tr(U t2 )) = Nf + (Tr(H t )) 

10 In contrast to the quenched case, where we do have flavour-singlet propagation at tree level. 


21 


















Cl HU)) 


(6.78) 


= ^Tr(t/ t )) — ~~~~ 


Let us define 


= ±±lnZ 
E ~ Nf dll " 1 


(6.79) 


for a theory with Nf quarks of equal mass. The partition function Z v is known explicitly in all 
generality 0. Up to an irrelevant overall factor the integration of eq. (|6.77|) gives 


N f 


Zv({ih}) = Iv+3-i(jh)]/ n^i “ Mj) > 

i>j 


(6.80) 


where I n (x) is the modified Bessel function, and the determinant in the numerator is over a matrix 
of size Nf x Nf (the indices i and j denote the matrix elements). In practice one may be mostly 
interested in the equal-mass case, where this expression simplifies considerably. Up to an irrelevant 
overall factor, 

Z v {{ii}) = det [I v+j -i((i)\ , (6.81) 

where again the determinant is taken over a matrix of size Nf x Nf. Then Ej,(/x) is known explicitly 
in terms of modified Bessel functions for all Nf. Next, using 




i 


2 N f \ 


(Ti (U + U^) 


(6.82) 


and the last of the identities ( 6.78 ), we obtain 


{T*m = 

NNi) = "/(^ + £) 


(6.83) 


Upon substituting eq. (3.83), the first two identities of eq. (6.75) as well as eq. (B.19) in eq. (6.74), 
we finally get 


(5°(.x)S°(0)) 
(P°(x)P°( 0)) 
(S a (x)S b ( 0)) 


_ qO (if) 


2 S 2 


s +-ur( iV /- 1 ) 


= cT ] + 


F 2 

2 E 2 

~F 2 


1M4 A (x) 

\Jb Zj 


■Nf 


l N f E v {fl) 1 EC(m) ( Ey(fi) t + 




N- 


f 


zab 


N 2 f~ 1 


C<y) s2 

U S ~-£p2 


2 Nf ^ + N f 




2 




A(x) 


E J 


-N f 


+ {2 Nj-F)- 2 + 


ir I 6N 2 — 4\ l T, v {ji) 


N f 


It 


A(x) 


(P a (x)P b ( 0)) = 


5 ab 


Nf ~ 1 


s 


4 F 2 


2 (Nj + 2) 


^U/f)\ 2 1 uC(m) 


E / 


AT 


+ 2JV|-8 + A£M4_5C 

fl 2-1 


A(x) 


(6.84) 


As a check of these expressions we can take the limit fi —> oo. For any finite z/, this freezes the zero¬ 
mode integral at U = 1, and we indeed find that the coefficients of the A(x) terms in the first three 
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equations of eq. ( [6.84 ) vanish in that limit. Similarly, the last coefficient approaches the required 
coefficient Y 2 /F 2 for reproducing the tree-level propagation of non-zero-momentum modes of the 
flavoured Goldstone bosons. 


The explicit expressions for the constant terms in eq. (6.84) are 

! ) 


r °(u) 

= s 2 ^ 

^O(v) 

Up 

= Z 2 N f 


Y 2 N f 

2 

4*° 

1 

II 


K~ loop 'M , Ar fzl- loop (M)V 

-v-+ N f -v 


1 Yl~ loo P{^) u 2 N f 


(¥)- 


s fi 2 

2 1 y}- 1o °p'(h) N f Yl~ loo P{fi) 


N 


f 




:i- i °op^) \ 2 


i 1 


(6.85) 


One notices that these QCD correlation functions at fixed topology need not be positive-definite. In 
the quenched theory the negative correlations had an immediate interpretation in terms of unitarity 
violation caused by eliminating the fermion determinant. Here, the origin of unitarity violation is only 
the restriction to fixed gauge field topology, and this peculiarity should disappear once we sum over 
v. Consider the main culprit, the constant part Cp of the (P°(x)P°( 0)) correlation function. For 
any given Nf we will have 

C° P {u) < 0 (6.86) 

for sufficiently large v. But when we sum over topology with the correct weight this does not occur. 
Consider the large-/r limit, where calculations can be made very explicitly. The distribution of winding 
numbers is then Gaussian ||], with average |33|] 

{v 2 ) = mYV/Nf = fi/N f . (6.87) 


In this large -n limit we can average over topology explicitly, keeping the leading order in 1 //j,: 


(C 0 p (u) ) 




Y 2 N f 

Y 2 N f 

0 . 


'l N f \ 

^ Y n 2 j 

1 ( v 2 )Nf 

/r [i 2 


( 6 . 88 ) 


To this order, negativity is precisely just avoided. This is consistent with the expression of eq. ( |6.76|) , 
where, after taking the limit U 1, we also find zero. 


As in the quenched theory, some of the above results have nice explanations in terms of general 
Ward identities, which are formally the same as in the quenched theory, with the obvious replacement 
N v —> Nf. Thus, from eq. ( 5.64 ), by combining again the identities corresponding to 0(x) = P°(x) 
and 0(x) = lo ( x ) we obtain to order e 2 , after integrating over space-time: 


C^y) 


Y 2 N f 


1 Yl~ lo °P{n ) 

/r Y 



(6.89) 
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which precisely matches what we found by explicit computations in eq. ( |6.85| ). Similarly, the identity 
for a non-singlet chiral flavour rotation of eq. (5.69|) requires 


CP = (AT|-l)±£Xi-^0u) 


(6.90) 


to 0(e 2 ) which again agrees precisely with the result of the explicit computation in eq. (6.85) 


Finally, we can also verify the Ward identity of eq. (5.64) for 0(x ) = P°(x ) directly. Using the last of 
the identities of eq. ( 6.7Sp we find, to this order in the e-expansion: 


^ (^{x)i^{x)) = - (Tr (U - t/ f )J) = 


(6.91) 


Then the Ward identity eq. (|5.64| ), after integrating over space-time, is seen to be identical to eq. ( 6.89 ), 
and hence also verified by our explicit calculations. 


There is no simple analogy of the quenched relation eq. (|5.71[), but the other relations of ref. [|30j] for 
Nf = 2 written in our normalization in eqs. ( 5.72[ ) and ( |5.73| ) should be satisfied. One first has to 
compute the chiral condensate in the Nf = 2 theory from eq. ( 6.79 ), and also find its first derivative. 
Next, plugging in Nf = 2 in eq. ( |6.85| ), and taking the p —» 0 limit, we find complete agreement. 
The cancellations required for this to happen are rather non-trivial, and hence constitute one more 
independent check on the calculation. 


7 Conclusions 

There are two very different regimes in which one can take the chiral limit of finite-volume QCD. The 
one that is normally considered is the one in which the volume is required to always exceed by far 
the Compton wavelength of all excitations, including the pseudo-Goldstone bosons that eventually 
become ultra-light. In this regime, there are still finite-volume effects, but they are exponentially 
small. This regime can be studied by lattice gauge theory simulations, but as the quark masses are 
decreased, it becomes increasingly difficult to provide the enormous volumes that are needed. This is 
where the usual p-expansion of chiral perturbation theory is relevant. For instance, if one wishes to 
measure the pion mass, the straightforward procedure is to ensure a posteriori that the linear extent 
of the lattice L by far exceeds the inverse pion mass 1 /m n , and then fit the Monte Carlo data to the 
usual zero-momentum projection of the pseudoscalar correlation function of eq. ( [5.46 ). Here we are 
advocating another procedure, which is much less costly in terms of computer resources close to the 
chiral limit. Instead of insisting on the large-volume condition L 1 /m n , one can go to any lattice 
size, and even go to the opposite regime L <C 1 /m n . In this regime another chiral perturbation theory, 
that of the e-expansion, is relevant. If one measures the same pseudoscalar correlation function in this 
regime, there is a modified formula that describes the correlation function. It depends on the same 
parameters of the chiral Lagrangian as in the infinite-volume case, and therefore it can just as well be 
used to determine these parameters. In particular, to lowest order it depends on the infinite-volume 
chiral condensate £, the pion decay constant F, and the quark mass m. By the Gell-Mann-Oakes- 
Renner relation, with corrections if required, this directly provides the pion mass m\ 7r . In this way the 
pion mass can be determined from a correlation function in a volume so small that not one single pion 
Compton wavelength can fit inside. 


In this paper we have discussed the new analytical formulae that arise in this regime of the e-expansion. 
In particular, we have derived the scalar and pseudoscalar correlation functions for quenched QCD in 
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sectors of fixed topological charge v, and for full QCD with Nf light flavours, also in sectors of fixed 
topology. Simultaneously we have discussed various aspects of the quenched e-expansion, especially 
its limitations, and the reason why it is particularly advantageous to consider the chiral limit of 
the quenched theory in sectors of fixed topology. We hope also to have made it apparent that the 
alternative replica formulation of quenched chiral perturbation theory can be very useful. 

As is well known, quenched QCD is a troubled theory, and in the end one should not consider our 
predictions for the quenched theory as more than a testing-ground for full QCD. In particular, we 
have found that the quenched finite-volume logarithms of the e-expansion do appear also in the 
scalar and pseudoscalar correlation functions. Here we have simply assumed that the quenched theory 
nevertheless can be defined over a certain range of scales. Our formulae can then be directly compared 
with results of quenched computer simulations. 

The predictions of Hansen [|3|] for full QCD, and our results here for full QCD in sectors of fixed 
topology, stand on a different footing. In full QCD the chiral e-expansion does not suffer from any of 
the pathologies of the quenched analogy. This is of course as expected, since in this case the theory is 
well-defined. Also here the analytical predictions can dramatically simplify certain aspects of Monte 
Carlo simulations, since they allow for a determination of the infinite-volume parameters of low-energy 
QCD from volumes much smaller than normally thought to be required. It will be interesting to see 
also these predictions of full QCD confronted with lattice gauge theory. 
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A Quenched Integrals over Zero-Momentum Modes 

In this appendix we give the technical details on how to derive all the G7(2|2) group integrals needed 
for the calculation of the flavoured correlation functions in section |5.l| . We have already outlined the 
idea in the main text. Here we also show how at this level one can switch between the supersymmetric 
and replica formulations. 


First a simple observation, most easily explained in the replica formulation. Since det J\A~ V Z v is a 
function of only, there exist a series of identities among zero-momentum mode expectation 

values. For example, all expectation values that can be obtained including j\4 sources such that 
both det A! and remain invariant, are identical. This gives us the first useful identity. Let 

X = EV(J\4 + J), with M. diagonal. Two sources introduced by either^ 


or 


J 


1 

EH 



J 


J_ ( 0 ij \ 

EH ^ 0 0 J ’ 


(A.l) 


(A.2) 


11 In the supersymmetric formalism the sources are added in the fermion-fermion block only. In the replica formalism 
they are added in the upper left-hand corner, and the matrix is then supplemented by N — N v = N — 2 additional 
mass-degenerate entries. 
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share the same det x and XX^ ■ This immediately leads to one useful identity, which is needed in section 

5H: 

(A-3) 


(ul + (c/-%) = 0 . 


Next we note that for a source matrix 


J = 


1 


ji 0 


(A.4) 


XV \ 0 32 ) 

we obtain the expectation value of eq. (|5.59|) after differentiating with respect to j\ and j 2 . What was 


evaluated explicitly in ref. 27] was rather the susceptibility 

5 2 


xM = 


SjiSj 2 


In Z L 


31=32=0 

/r 2 (/W(/x ) 2 - K u+ 1 (v)K u - i( M )) ~ Iu+iWv-l 


(A.5) 


We therefore have, after substituting from eq. ( 4.21 ): 

<5 2 


Sj\Sj 2 


ji=ja=0 


\((U 11 + U^ 1 ){U 22 + u£)) 


= Xv(fJ-) + 


V 

1 H-2 J 




(A.6) 


where in the last line the remarkable simplification is due to a series of Bessel function identities. 


To see how we can use this result, and that of eq. ( 5.38 ), to evaluate all needed expectation values, 
let us, as an example, consider a source matrix of the form 


J = 


This provides us with the expectation value 


J_ ( 0 3 

El/ l 0 0 


((C/21 + (CC-%) 2 } = 4 


sr 


(A.7) 


(A-8) 


3=0 


To evaluate it, we note that the square roots of the eigenvalues of the corresponding matrix XX' with 
X = M. + J are 


mr = -^(2g 2 + J 2 + (/ + 4m 2 J 2 ) 1 / 2 ) 1/2 
m 2 = -j= (2/j? + j 2 - (j 4 + 4m 2 j 2 ) 1/2 ) / 


(A.9) 


and, in the replica formalism, N — N v = N — 2 additional degenerate eigenvalues fi. By the chain rule, 
and using that (detx)~ u Z u clearly is a symmetric function of the eigenvalues mi j2 , we get 


dj 2 


1 d 


3 =0 


2 n dmi 


1 92 Z 

mi=m-2=[i ^ dm 1 


1 


d 2 


777-1 —777-2 —/-C 


2 dm\dm 2 


Z h 


(A.10) 


777-1 —7712—/i 
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where on the right-hand side Z v is the partition function of diagonal sources 


M + J = 


1 j mi 0 

Sh 0 m 2 


(A.ll) 


Performing the required differentiations, and substituting eq. ( 5.38 ) and eq. ( 5.4C| ), we find the compact 
expression 


((t/21 + ( U - 1 )l 2 ) 2 } = 2 


1 

II S E 


(A.12) 


Next, consider a purely imaginary source matrix such as 

1 ( iji 0 


J = 


El/ l 0 ij 2 


(A.13) 


With x = E1/(A4 + J), the eigenvalues of XX^ are either fi 2 + j 2 . fj 2 + j 2 (in the supersymmetric 
formulation), or /x 2 + . 7 2 ,/x 2 + j 2 ,{^ 2 } (hr the replica method). In the latter formulation there are 
N — N v degenerate eigenvalues in the last set. In both cases we have 


det(x) = 1 + 


in the quenched limit. Now consider 

d 2 


<9ji <9.72 


— [(det x) v Z„\ h=j2=Q = - — - 


Ui 


v 2 2iv d 


1 + 


m 




l L dji 


Z„ 


d 2 „ 

+ o • o • Z h 

31 = 32=0 


(A.14) 


(A.15) 


J1=J2=0 


Here the right-hand side actually vanishes, as follows from switching to the square roots of the eigen¬ 
values of xx^- We thus have 


d 2 


djidj 2 


Z, 


v 2 2 iv d _ 

= “W H- TT^Z, 


j 1=J2=0 




M <9.71 


(A.16) 


jl=j2=0 


To evaluate the remaining single derivative, we can proceed analogously, using 


o = |_[(detx)-^W „ = 


jl=j2=0 


(A.17) 


Putting everything together, we finally have 


<(c/n - (T/- 1 ) 11 )(C/ 22 - (t/- 1 )^)} = - 4 ^- 


<9j2 




4zx 2 

/^ 2 


j 1=J2=0 


(A.18) 


These, and those given in section |5|, are all the expectation values we need for evaluating the flavoured 
correlation functions. 
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B U(N ) Group Integral Identities 


In section |we used some identities for group integrals over U(Nf) of the zero-mode action eq. (6.77). 
These identities, and others needed if one wishes to go beyond leading order, can be viewed as 
Schwinger-Dyson equations on that particular group manifold. Here we outline how to derive such 
identities. 


To begin, we introduce the notion of left- and right-handed differentiation on the group. Let t a be the 
generators of the algebra u(Nf) in a given representation. Differentiation on group elements of U(Nf) 
in the same representation can be defined by either 

F(Ue ie “ ta ) = F{U) + e a V a R F(U) + ... , (B.l) 

for right-handed differentiation, or 

F(e ie “ ta U) = F(U) + e a V a L F(U) + ... , (B.2) 

for left-handed differentiation. The latter is probably the most intuitive to work with, and we will use 
it in what follows. For simplicity of notation we omit the subscript L from now on. One immediately 
sees that an explicit representation of V a is 

V “ = ■ (B ' 3) 

By hermitian conjugation, or by invoking unitarity, one finds the way V a acts on hermitian conjugates. 
For example, 

V a L rt = - iUh a . (B.4) 

The derivative in eq. ( |B.3| ) satisfies the Leibniz rule 

V a (FG) = (V a F)G + F{V a G) , (B.5) 


and the Lie algebra 

[V a ,V fe ] = / a6c V c , 

where f abc are the structure constants. 

Left-invariance of the Haar measure on \J(Nj) implies that 

JdU V a F{U) = 0 , 


(B.6) 


(B.7) 


which, in conjunction with the Leibniz rule, leads to a simple rule for partial integration. 


Let us now focus on the zero-mode theory of (6.77). Schwinger-Dyson equations of this theory are 
obtained by applying the rule (B/i) through an insertion of a set of functions of U in the group integral. 
As an example, consider 


0 = 


dU Ti[t a \7 a (F(U)(detU) u e s W)\ 


(B.8) 


where 

S(U) = 'r(U + tf) . 


(B.9) 
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We normalize the generators by Tr (t a t b ) = \5 ah , and we also need the U(Nf) completeness relation 
= \8ii8jk in that normalization. Choosing F{JJ) = U in ( |B.8 ) leads to the identity 


(Tr(17 2 )) = N f - — (Tr(17)) , (B.10) 

which was used repeatedly in section |6|. Because the integration measure in ( p.8| ) is not real when 
i//0, the analogous identity for the hermitian conjugate does not follow trivially from this. Rather, 
from choosing F(U ) = U' in ( p.8| ) one finds 

(Tr((b rt ) 2 )) = N f + 2{U ~^ Nf) (Tr(f/t)) . (B.ll) 

As a consistency check we note that it can also be obtained from eq. ( p.lOj ) by taking the hermitian 
conjugate while simultaneously letting v —> —v. Next, by considering an identity such as 

0 = JdU V°[Tr {U - 17 t )(det U) u e S{JJ) ] , (B.12) 

one finds 

|([Tr(C/-C/ t )] 2 ) = - (Tr{U + £/*)) - vN f (Ti(U - t^)) . (B.13) 

In order to simplify the last term, one can consider the identity obtained from 

0= jdUV°[(detUyeW], (B.14) 


which gives 


Tr (U-tf)) = 

/ fi 


Inserting this into eq. (B.13) gives 

|([TV(t/-t/t)] 2 ) = -2 + 


(B.15) 


(B.16) 


Let Z v be as defined in eq. (5.77). Then the relation 


with Sj,(/i) as defined in eq. (|6.7S ). We also need to extract individual terms on the left-hand side. 

(B.17) 


i([TV(t/ + t/t)] 2 } = 


Z v <9/x 2 


gives 


([Tr(T/ + [/t)] 2 ) = 4JV, 


(B.18) 


Adding and subtracting ( p.lG| ) and ( p.l^ ) hnally provides the needed identities: 

'KM 


Tr(L r )Tr(t/ t ) 
(Tr(t/)) 2 + (Tr([/t)) 2 


= Nt 


= 2 N f 


E 

KM 


+ n JM^Y + 1 -EM-Z1 

V s ) II E fi 2 

+ A RML) 2 + L^iM4 

V E J jjL- n E 


(B.19) 
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